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An ab initio relativistic coupled-cluster theory of dipole and quadrupole 
polarizabilities: Applications to a few alkali atoms and alkaline earth ions 
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We present a general approach within the relativistic coupled-cluster theory framework to calculate 
exactly the first order wave functions due to any rank perturbation operators. Using this method, 
we calculate the static dipole and quadrupole polarizabilities in some alkali atoms and alkaline earth- 
metal ions. This may be a good test of the present theory for different rank and parity interaction 
operators. This shows a wide range of applications including precise calculations of both parity 
and CP violating amplitudes due to rank zero and rank one weak interaction Hamiltonians. We 
also give contributions from correlation effects and discuss them in terms of lower order many-body 
perturbation theory. 

PACS numbers: 31.15.Ar,31.15.Dv,31.25.Jf,32.10.Dk 



For a long time, studies of dipole and quadrupole po- 
larizabilities have been highly interesting in many im- 
portant aspects for both neutral atoms and ions [l], Q, d, 
0, % @, 0, S, ©, EH El ■ These quantities are used in the 
case of ions to remove the quadratic Stark shifts in the re- 
cent proposed optical frequency standards where suitable 
instruments are not available for very precise measure- 
13| . Relativistic effects are also investigated 



12, 



ments 

for these properties Elj]. Often theoretical studies 
are limited by many approximations to only approximate 
accuracy of the results 17], |2l|. In many cases, semicm- 
pirical methods, which combine the calculated El and 
E2 matrix elements of important states and the exper- 
imental excitation energies, are used to calculate these 
quantities [18|, Ell- Core electron contributions are al- 
ways only estimated in this approach [l], EH- Therefore, 
an ab initio method is necessary to be able to, in prin- 
ciple, explain the importance of the electron correlation 
effects in the above properties and also test the many- 
body theories employed. These results are largely used 
to determine the van der Waals coefficients or dispersion 
factors 0,0, EH . l22l |. In some cases, the Dirac-Fock 
(DF) method and relativistic random-phase approxima- 

231 or scalar 



tions (RPA) for the Dalgarno-Lewis [2l|, [22 
relativistic Douglas-Kroll [TBI, EB| effective Hamiltonians 
are used to determine these quantities. A number of cal- 
culations on these quantities using molecular codes and 
pseudopotentials have also been reported [HEM 3 11 \- 



Using the relativistic coupled-cluster (RCC) theory, we 
report here a novel approach for the first time to calculate 
wave functions up to all orders in the residual Coulomb 
interaction and first order due to any kind of perturbed 
operator in one valence systems. This takes care of the 
sum-over-states approach of the many-body perturbation 



theory (MBPT) as an exact solution. This method can 
be applied to calculate the first order wave function due 
to dipole and quadrupole transition operators. Hence, 
polarizabilities can be obtained by calculating the expec- 
tation values of the corresponding operators. As an ab 
initio test of the theory, it is also possible to calculate the 
light shift ratio in ions like Ba + using the present method 
as proposed in [24| • Moreover, parity non-conserving 
(PNC) and CP violating amplitudes can be calculated 
very precisely using this method due to rank zero or rank 
one weakly interacting Hamiltonian, which are the most 
challenging studies during the last three decades [H, W 



To demonstrate the application of the method, we con- 
sider three systems each from alkali atoms (Li, Na and 
K) and alkaline earth ions (Be + , Mg + and Ca + ) to cal- 
culate their dipole and quadrupole polarizabilities having 
different angular momentum selection rules. 

The energy shift, AE(J n , M n ) of any state |J„M„ >, 
with principal quantum number n, in a direct current 
(dc) electric field E = £z can be expressed as 

AE{J n ,M n ) = ~a\j n ,M n )£ 2 , (0.1) 

where a 1 (J„,M„) is defined as the static polarizability 
of state \ J n M n >. Further, a 1 (J„,M„) can be divided 

as 

3M2 - J„(J„ + 1) 



a 1 (J„,M„) = al(J n ,M n ) + 



-a L 2 (J n ,M n (J).2) 



J n (2J n - 1) 

Here a^J^, M n ) and a\{J n , M n ) are known as the scalar 
and tensor polarizabilities, respectively. From the first 
order perturbation equations, these parameters can be 
expressed as the sum over intermediate states 

a 4 ( J n , M n ) = -2 2^ , (0.3) 



E n — Ef. 
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where i represents either or 2, C\ are the appropriate 
angular co-efficients, z is the z component of the position 



2 



vector r and E's are the unperturbed energy levels. Since 
z can be expressed in terms of the spherical harmonics of 
rank one (Y 10 (9,4>)), the above matrix elements will be 
non-zero between opposite parity states satisfying the El 
transition selection rules. The o^'s can be expressed in 
terms of the reduced matrix elements of the El operator 
(D = er) as follows 



3(2J„ + 1) 



k^n 



E„ — Ei, 



(0.4) 



where |$„) is the open-shell reference state con- 
structed by appending the valence electron 'n' orbital to 
the closed-shell Dirac-Hartree-Fock (DF) wave function 
(|$o); i- e - l^n) = a nl^o))- I n the above expression, the 
coupled-cluster (CC) excitation operators, T excites only 
the core electrons and S n excites either only the valence 
electron 'n' or along with necessary core electrons. 

To get both the unperturbed and perturbed wave func- 
tions of Eqn. (|0.8|> . we express the T and S n operators 
by 



and 



a\{J n ) 



40J(2J-1) 



3(2J + 3)(2J+ 1)(J+ 1) 



1/2 



J 71 



1 Jk 
J 2 



y^(-i) ,7 " +j '° +i 

k^n 
\(Jn\\D\\J k )\* 



T = T (0) + T (1) 

Sri — sL ~k sL ■ 



(0.11) 
(0.12) 



1 Jn 2 J E n — Ek 
Similarly, the static quadrupole polarizability can be 



(0. 



expressed as 

a 2 {J n ,M n ) 



k^n 

-2 



5(2J„ 



\(J n M n \Q\J k M k )\ 2 

E n — E k 
\- |(Jn||Q||J fc )| 2 

E„ — Ek 



1) ^ 

1 k^n 



, (0.6) 



where Cq is the corresponding angular factor and Q = 
— |(3z 2 — r 2 ) is the E2 operator which has different se- 
lection rules than the El operator. 

The above expressions for both the polarizabilities can 
be expressed in a general form as 

a(J n M n ) = <* n |0|*„> 

= (*M\0\*W) + (*W\D\*W) 

= 2<*W|O|*(0>>, (0.7) 

where the exact wave function of the nth state can be 
written in terms of the original atomic wave functions and 
first order corrections due to the corresponding dipole or 
quadrupole operators 0{= D or Q); i.e. 

|¥ n ) = + (o.8) 

Note that the C angular factors from the a's are absorbed 
in the first order wave functions. 

We show that it is possible to calculate a(J n M n )'s 
exactly by calculating both the \^n^) and j*! ) using 
single many-body theory which can avoid the sum-over- 
states approach given above. In our approach, we would 
like to obtain the first order perturbed wave function as 
a solution to the following equation 

(H( DC } - E^)\9^) = - Hg)|*8»)>, (0.9) 

where Hq 00 ' and H; nt are the Dirac-Coulomb (DC) and 
interaction Hamiltonians due to El or E2 operators, re- 
spectively. The En and En are the zeroth and first 
order energies of the nth state, respectively. 

In the RCC theory, we express the exact wave function 
for one valence (n) state of a system as 

|*n> = e T {l + S n }\<S> n ), (0.10) 



where and S n 0] are the CC operators for the DC 
^Hamiltonian andTW and are the corresponding first 
order excitation operators due to the interaction Hamil- 
tonian. To calculate aj(J„M„), only linear terms involv- 
ing T^ 1 ) or Sn^ operators are kept from the exponential 
function of Eqn. (|CTT0|) 



|*n) = e Ti0) {l 



(l+SW)TM+SM}\* n )(0.13) 



Hence, the unperturbed and perturbed wave functions 



can be separated as 

|*(°)) = ^ (0, {i + <#)}|$„), 



(0.14) 



'{(l + SW)TW + sW}|* n ). (0.15) 



We consider only the single and double excitations 
from the RCC method (CCSD method) in our calcula- 
tions, which is proved successful to cope with the electron 
correlation effects in most one valence systems; i.e. 



T = T 1 +T 2 

S n = S\n ~\~ S 2 n- 



(0.16) 
(0.17) 



First we solve the unperturbed and Sv ampli- 
tudes by solving usual CC equations, then these ampli- 
tudes are used to determine the and Sv amplitudes 
in the following equations 



<$5|f4 DC) T(l)|$ ) 



($5|H int |$ ) (0.18) 



^|(H^-IP„)5«|$„) = 



-(<S> n \ [H int {l + S<°>} + H^ C) T«] |$„), (0.19) 

where the subscript N represents normal order form of 
the operators, IP„ is the ionization potential of the nth 
state and the symbol * represents excited states with 
respect to the corresponding reference states. 

After getting both the unperturbed and perturbed am- 
plitudes, we evaluate polarizabilities using the following 
expression 



3 



ai(J n M n ) 



{$ n \{l + S n }e T 'Oe T {l + S n }\<t> n ) 
(<S> n \{l + Si}e Ti e T {l + S n }\$ n ) 

< 9> n \S^0^{l + S n 0) } + {1 + S* 0)t }OC°fc« + {1 + 4 0)t }(T« t oW + 0C°)T«){1 + S n 0) }\^ n 



(6.20) 



TABLE I: Static dipole and quadrupole polarizabilities in al- 
kali atoms: Li, Na and K. 



TABLE II: Static dipole and quadrupole polarizabilities in 
alkaline earth-metal ions: Be+, Mg+ and Ca+. 



Atoms Expts 



Others 



This work 



Atoms Expts 



Others 



This work 



Dipole polarizabilities 



Dipole polarizabilities 



Li 



Na 



K 



Li 



Na 



K 



164(3.4) [28] 
148(13) [29] 



163.73 [15], 164.6 [9J 
165.01 [17], 164.111 [6] 



162.29 



162.; 



286.01 



159.2(3.4) [28] 163.07 [18], 164.89 [15] 
164.6(11.5) [29] 160.7 [9], 165.88 [17] 

292.8(6.1) [28] 290.1 [18], 289.5 [9] 
305(21.6) [29] 301.28 [15], 285.23 [17] 

Quadrupole polarizabilities 

1424(4) [21], 1393 [9] 1421.28 
1423.266(5) [6], 1424 [30] 
1423 [31] 

1885(26) [21], 1796 [9] 1899.67 
1878 [30], 1879 [31] 

5000(45) [21], 4703 [9] 4919.71 
5000 [30], 5001 [31] 



Be 4 



Me 



Ca 4 



Be 4 



Ms 



Ca 4 



34.62(26) [4] 
33.0(5) [35] 
33.8(8) [36J 

70.89(15) [4] 
75.3(4) [35] 



24.93 [9], 24.63 [32] 
25.04 [33], 16.74 [34] 

34.144 [4], 38.7 [37] 
33.68 [9], 37.2 [34] 

34.0 [38], 38.9 [39] 

70.872 [4], 112.4 [3JJ 

71.01 [9], 87(2) [3] 



72.5(19) [35] 76.9 [38], 96.2 [34] 

Quadrupole polarizabilities 

52.93 [9], 55.42 [33] 
55.71 [34], 52.4 [2] 

150.2 [9], 187.66 [34] 
150.15 [38] 

1171 [9], 727.55 [34J 
1303.51 [381 



24.11 



34.59 



73.86 



53.80 



156.17 



706.59 



where for computational simplicity we define CK°) = 



Oe 1 



and N n 0) 



On & 



S n . We com- 



pute 0(°) in two steps as effective one-body and two-body 
terms and substitute in the above equation. We account 
for contributions from the normalization factor expressed 
as 



Norm. = <*„|0|tf„){ 



1 + JVn 



!}• 



(0.21) 



The above approach can easily be extended for a very 
precise calculations of PNC amplitudes [25} and CP vi- 
olating electric dipole moments (EDMs) in atoms [26j 
by obtaining the first order wavefunctions due to nuclear 
spin independent and dependent weak interaction Hamil- 
tonians in the place of dipole and quadrupole operators. 

We construct the relativistic single particle orbitals us- 
ing Gaussian type orbitals (GTOs) and we consider the 
finite size of the nucleus assuming a Fermi charge distri- 
bution as discussed in 27J . One can use length and veloc- 



ity gauge expressions for El and E2 operators to verify 



accuracy of the results. However, we have used only the 
length gauge expressions which are more stable than oth- 
ers in our calculations. We present our results in Tables 
U and UH and compare them with the experimental re- 
sults and other calculated values. Numerous calculations 
are available for the dipole polarizabilities in the neutral 
systems. However, we present only a few recent calcula- 
tions in these tables. To our knowledge, there are only 
a few experimental results for the dipole polarizabilities 
whereas no results are found for the quadrupole polariz- 
abilities, though there exist a number of calculations. 

We present the DF and different RCC results in Table 
IIIII to understand the role of electron correlation effects. 
As seen from this table, OS[ n ^ and its complex conju- 
gate (cc) terms contribute predominantly. To interpret 
these contributions, we break down this RCC term into 
some of the leading order MBPT diagrams as shown in 
the Fig. ma)). It is obvious from this relationship that 
the above term consists of dominant DF (Figfjjai)), core- 
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TABLE III: Contributions from DF and important perturbed 
CC terms for the dipole and quadrupole polarizabilities. 



Atoms 


DF 


OS^+CC 


OS^+cc 


Norm. 


Others 




Dipole polarizabilities 






Ll 


168.95 


164.21 


-0.40 


n no 
— 0.08 


1 A A 

— 1.44 


JN a 


188.17 


169.09 


-1.26 


n on 
— 0.29 


—4.00 


T/ 

rv 


398.15 


313.32 


-6.93 


O 1 1 

— Z.ll 


1 o oo 
— 18.28 


Be 


24.81 


24.33 


-0.08 


—0.01 


—0.13 




38.39 


35.57 


-0.43 


—0.04 


—0.51 


Ca + 


94.62 


79.84 


-1.48 


—0.32 


—4.18 




Quadrupole polarizabilities 






Li 


1484.98 


1444.87 


0.00 


-0.72 


-22.87 


Na 


2230.62 


2012.95 


0.00 


-3.34 


-109.94 


K 


7099.70 


5562.48 


0.00 


-36.38 


-606.39 


Be+ 


54.97 


54.23 


0.00 


-0.02 


-0.41 


Mg+ 


171.34 


161.09 


0.00 


-0.18 


-4.74 


Ca+ 


952.65 


750.39 


0.00 


-3.07 


-40.73 



particle physics. We also investigated electron correlation 
effects from the dipole and quadrupole polarizabilities 
calculations in six different systems and highlight their 
behavior. We observed that higher order correlation ef- 
fects become important in large systems; suggesting that 
the method can be employed rigorously in heavy systems 
to obtain accurate results, and can be studied as ab initio 
tests. Indeed, this approach can also be extended to de- 
termine frequency dependent (dynamic) polarizabilities 
which we defer to our next studies. Conclusively, the 
many-body aspects presented here would certainly be of 
interesting for both physicists and quantum chemists. 

We are grateful to Professor B. P. Das for many useful 
discussions and suggesting us to carry out this work. We 
are also delighted to thank Professor P. Fulde for very 
enlightening discussions and hospitalities at MPI-PKS. 




FIG. 1: Break down of important perturbed CC diagrams into 
some of the important lowest order MBPT diagrams. The V 
represents Coulomb operator and a and p represent, occupied 
and unoccupied orbitals, respectively. 



polarization (FigQJaii) and [TJaiii)) and pair-correlation 
(FigOJaiv)) effects, therefore this results in large contri- 
butions for all the systems. However, we also show in the 
Fig. |]Jb)) another important RCC term (OS!^) as the 
sum of different types of core-polarization effects. From 
Table Hill we find that contributions from other higher or- 
der terms are non-negligible while going towards the large 
systems. We also notice from this table that the trend of 
the correlation effects in neutral and ionic systems with 
the same electronic configurations are different. 

In conclusion, we have developed for the first time a 
general method in the relativistic coupled-cluster theory 
to calculate the first order wave functions due to any 
rank operators for both the parities. This technique can 
be applied to diverse areas of physics ranging from polar- 
izabilities to probe physics beyond the Standard Model of 
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